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1.  INTRODUCTION 

In  late  April  and  early  May  1986,  Air  Force  Geophysics  Laboratory  (AFGL) 
personnel  from  the  Atmospheric  Optics  Branch  participated  in  a  coordinated  field 
program  at  Pennsylvania  State  University  (PSU)^.  The  purpose  of  this  program  was,  in 
part,  to  compare  altitude  profiles  of  Cn,  the  refractive  index  structure  parameter^, 
obtained  from  different  measuring  devices  and  to  use  the  acquired  data  for  checking  the 
reliability  of  certain  models  for  Cn-  Profiles  of  this  parameter  can  be  obtained  using  a 
stellar  scintillometer.  The  one  used  at  PSU  by  AFGL  is  a  modification  of  that  originally 
developed  by  Ochs  et  aF.  This  instrument  measures  the  variance  of  stellar  intensity  for 
a  first  magnitude  or  greater  star.  This  variance,  by  use  of  spatial  filters,  can  be 
converted  into  a  spatially  and  temporally  averaged  profile  of  Cn  for  7  slant  path  levels 
where,  by  slant  path  level,  we  mean  distance  measured  along  the  line  of  sight  from  the 


(Received  for  publication  10  August  1990) 

^Markson,  R.  M.,  Anderson,  B.  W.,  Fairall,  C.  W.  ,  Thomson,  D.  W.,  White,  A.  B. 
and  Syrett,  W.  J.  (1989),  Atmospheric  Turbulence  Measurements  in  Support  of  Adaptive 
Optics  Technology,  RADC-TR-89-289,  ADA220445. 

^Hufnagel,  R.  E.  (1978),  Propagation  Through  Atmospheric  Turbulence,  in  The 
Infrared  Handbook,  USGPO,  Washington  D.  C.,  Chap.  6,  1-56. 

^Ochs,  G.  R.,  Wang,  T.,  Lawrence,  R.  S.  and  Clifford,  W.  S.  (1976),  Refractive 
Turbulence  Profiles  Measured  by  One-Dimensional  Spatial  Filtering,  Appl.  Opt.,  15:  2504-2510. 
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instrument  to  the  chosen  star.  Because  the  elevation  angle  for  the  star  used  is  known  as 
a  function  of  time,  the  slant  path  distances  can  be  converted  to  obtain  the 
corresponding  altitude  above  ground. 

The  Cn  values  obtained  for  adjacent  levels  are  not  entirely  instrument 
ind('p('ndent  because  there  will  be  some  overlap  in  the  weighting  functions  used  to 
separate  the  spatial  frequencies.  It  is  usually  assumed  that  readings  from  levels  1,  4,  and 
7  are  instrument  independent.  Figure  1  shows  a  typical  scintillometer-derived  profile'  for 
Cn-  The  scale  shown  has  been  chosen  to  show  the  usual  altitude  range  of  interest  for 
this  parameter  and  its  typical  variation.  We  see  that  the  scintillometer  does  not  provide 
information  regarding  Cn  for  high  and  low  altitudes  but  does  give  a  profile  from  about  2 
km  to  17  km.  It  should  be  emphasized  that  the  “points”  on  this  graph  really  repr('sent 
values  centered  at  the  indicated  altitude  but  come  from  averaging  over  distances  of  the 
order  of  a  kilometer  and  a  time  duration  of  2-3  minutes. 

It  has  been  brought  to  our  attention  by  one  of  the  co-investigators  in  this 
measurement  program  from  Rome  Air  Development  Center  (R.4DC)''  that  correlation 
coefficients  for  Cn  values  derived  from  AFGL  data  (192  profiles)  between  different 
levels  seem  surprisingly  large.  Furthermore,  this  phenomenon  has  been  observed  for 
other  scintillometer  data  sets  from  other  measurement  programs  and  the  conclusion  is 
that  the  assumption  that  the  values  of  Cn  from  levels  1,  4,  and  7  are  independent  is 
incorrect.  RADC  also  used  a  stellar  scintillometer  at  PSU  but  their  data  (140  profiles) 
does  not  lead  to  such  a  high  level  of  correlation. 

In  the  absence  of  any  other  considerations,  the  correlation  coefficients  calculated 
by  RADC  for  the  AFGL  scintillometer  data  would  indicate  that,  for  all  levels,  Cn 
values  are  instrument-dependent,  contradicting  the  above  stated  assumption  of 
independence  between  levels  1,  4,  and  7.  It  is  the  purpose  of  this  report  to  further 
consider  these  observations  and  to  suggest  other  possible  reasons  for  the  high  degree  of 
correlation  observed. 

2.  CALCULATIONS 

Data  is  available  for  the  evenings  of  30  April,  2  May,  3  May,  4  May,  5  May  and 
6  May.  The  number  of  profiles  for  each  evening  are  12,  39,  5C,  2C,  16,  and  43 
respectively,  a  total  of  192. 

“'Stebbins,  D.  W.  (1987).  Analysts  of  the  Inter-layer  Independence  of  Stellar 
Scintillometer  Profiles  of  Cl,  FtADC-TH-87-l.b6.  ADA189372. 
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2.1  Altitude  versus  Level  Number 

We  began  by  looking  at  the  difference  between  slant  path  levels  and  altitude 
above  ground  for  the  AFGL  scintillometer  profiles.  During  the  PSU  program,  the 
instrument  was  pointed  at  the  same  star  (Arcturus)  at  a  near  zenith  position.  The 
actual  altitude  for  each  measurement  was  calculated.  Then  the  mean  and  standard 
deviation  for  all  measurements  were  calculated.  (See  Table  1.)  Because  the  standard 
deviation  is  of  the  order  of  5  percent  of  the  mean  for  each  slant  path  level,  we  shall 
assume  that  the  Cn  values  from  a  given  level  correspond  to  the  same  altitude. 

2.2  Correlation  and  Regression 

The  methods  of  linear  regression  and  correlation^  are  frequently  used  to 
inv(^stigate  the  relationship  between  two  (or  more)  variables.  By  the  method  of  least 
squares  the  best  straight  line  fit  between  two  data  sets  can  be  obtained.  The  linear 
correlation  coefficient,  r,  measures  the  strength  and  direction  of  the  linear  association  of 
the  two  variables.  The  square  of  the  correlation  coefficient  is  the  fraction  of  the 
variation  of  one  variable  that  is  explained  by  least  squares  regression  on  tin'  other 
variable.  It  should  be  noted  that  correlations  based  on  averages  are  usually  too  high 
when  applied  to  individual  values®. 

We  have  calculated  correlation  coefficients  between  each  level  for  both  Cn  tind 
log(Cn)  for  each  evening  and  for  the  entire  program.  These  results  are  presented  in 
Tables  2-8.  Each  table  presents,  in  matrix  form,  two  sets  of  correlation  coefficients: 
rji(i,j)  and  r£^(i,j)  for  the  correlations  of  the  numerical  and  logarithmic  values  of  Cn 
between  level  i  and  j.  Tables  2-7  are  for  each  evening  and  Table  8  contains  resvdts  for 
the  entire  measurement  program. 

When  we  look  at  the  results  for  all  192  profiles,  we  see  for  both  Cn  and  log(Cn)  a 
significant  positive  correlation  between  all  levels.  Correlation  coefficients  decrease  as 
the  altitude  separation  between  levels  increases.  Level  1  (at  2  km)  and  level  2  (at  3.1 
krn)  are  the  closest  and  show  the  largest  correlation.  Levels  1  and  7  (at  16.6  km)  show 
the  least. 

For  the  entire  data  set  a  picture  much  different  than  that  of  any  individual 

®Sce,  for  example,  Brockett,  P.  and  Levine,  A.  (1984),  Statistics  and  Probability  and 
Thf'ir  Applications,  chapter  9,  Saunders  College  Publishing. 

®For  a  detailed  example  of  this  phenomenon  see  Moore,  D.  and  McCabe.  C.  (1989). 
Introduction  to  the  [Practice  of  Statistics,  p.  210,  W.  H.  Freeman. 
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Table  1.  Mean  altitude  and  standard  deviation  for  each  level. 


Level 

Mean  Altitude  (Km) 

Standard  Deviation  (Km) 

1 

1.97 

0.11 

2 

3.05 

0.17 

3 

4.66 

0.26 

4 

6.54 

0.36 

5 

8.42 

0.46 

6 

12.54 

0.69 

7 

16.57 

0.91 

Table  2.  Correlation  coefficients  between  levels  for  Cn:  30  April  1986  (12  profiles). 


Correlation  coefficients  between  levels  for  log(Cn):  30  April  1986  (12  profiles). 


1 

2 

3 

4 

5 

6 

7 

1 

0.96 

0.95 

0.72 

0.70 

0.18 

2 

0.99 

0.57 

0.63 

0.15 

3 

0.61 

0.66 

0.22 

4 

0.86 

0.09 

0.58 

5 

-0.04 

0.62 

6 

0.02 

7 

5 


Table  3.  Correlation  coe  fficients  between  levels  for  Cn:  2  May  1986  (39  profiles). 


1 

2 

3 

4 

5 

6 

^  7 

1 

0.86 

0.43 

-0.26 

-0.25 

0.17 

0.08 

2 

0.71 

-0.17 

-0.18 

0.26 

0.12 

i  3 

0.43 

0.28 

0.34 

0.06 

4 

0.88 

0.05 

0.14 

5 

-0.06 

0.24 

6 

0.24 

7 

Correlation  coefficients  between  levels  for  log(Cn);  2  May  1986  (39  profiles). 


1 

2 

3 

4 

5 

6 

7 

:  1 

0.81 

0.44 

-0.31 

-0.28 

0.04 

0.09 

i  2 

0.72 

-0.31 

-0.31 

0.15 

0.10 

3 

0.22 

0.08 

0.29 

0.04 

4 

0.90 

-0.12 

0.05 

5 

-0.15 

0.15 

6 

0.23 

7 

I  j 

6 


Table  4. Correlation  coefficients  between  levels  for  Cn:  3  May  1986  (56  profiles). 


1 

2 

3 

4 

5 

6 

7 

1 

0.92 

0.85 

0.39 

0.30 

0.24 

0.27 

2 

0.93 

0.42 

0.25 

0.46 

0.38 

3 

0.65 

0.42 

0.31 

0.34 

4 

0.87 

0.10 

0.36 

5 

0.12 

0.34 

6 

0.60 

7 

Correlation  coefficients  between  levels  for  log(Cn):  3  May  1986  (56  profiles). 


5 

6 


0.22 


Table  5. Correlation  coefficients  between  levels  for  Cn:  4  May  1986  (24  profiles). 


1 

2 

3 

4 

5 

6 

7 

1 

0.33 

-0.03 

-0.07 

0.20 

0.29 

0.18 

2 

0.73 

0.34 

0.06 

-0.14 

-0.55 

3 

0.65 

0.11 

-0.24 

-0.54 

4 

0.41 

-0.03 

-0.15 

5 

0.35 

0.17 

6 

0.37 

7 

Correlation  coefficients  between  levels  for  log(Cn):  4  May  1986  (24  profiles). 


1 

2 

3 

4 

5 

6 

7 

1 

0.19 

-0.06 

0.01 

0.21 

0.17 

-0.05 

2 

0.57 

0.39 

-0.04 

-0.18 

-0.49 

3 

0.81 

0.21 

-0.19 

-0.44 

4 

0.45 

0.05 

-0.13 

5 

0.23 

-0.03 

6 

0.17 

_ 

_ 1 

Table  7. Correlation  coefficients  between  levels  for  Cn^  6  May  19S6  (43  profiles). 


1 

2 

3 

4 

5 

6 

7 

1 

0.88 

0.85 

0.80 

0.83 

0.50 

0.25 

2 

0.99 

0.86 

0.64 

0.28 

0.27 

3 

0.89 

0.62 

0.28 

4 

0.71 

0.40 

0.27 

5 

0.63 

0.12 

6 

0.04 

7 

Correlation  coefficients  between  levels  for  log(Cn):  6  May  19SC  (43  profiles). 


1 

2 

3 

4 

5 

6 

7 

1 

0.77 

0.72 

0.52 

0.66 

0.66 

0.17 

2 

0.88 

0.73 

0.62 

0.43 

0.21 

3 

0.88 

0.68 

0.38 

0.20 

4 

0.69 

0.36 

0.13 

5 

0.57 

0.04 

6 

0.08 

7 
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Table  8.  Correlation  coefficients  between  levels  for  Cn:  30  April-6  May  1986  (192  i)rofiles). 


1 

2 

3 

4 

5 

6 

7 

1 

0.91 

0.82 

0.44 

042 

0.40 

0.37 

2 

0.93 

0.52 

0.52 

0.43 

3 

0.72 

0.44 

0.40 

4 

0.87 

0.34 

0.46 

5 

0.42 

0.56 

6 

0.60 

7 

Correlation  coefficients  between  levels  for  log{Cn);  30  April-6  May  1986  (192  profiles). 


1 

2 

3 

4 

5 

6 

7 

1 

0.88 

0.76 

0.54 

0.61 

0.62 

0.52 

2 

0.92 

0.66 

0.64 

0.43 

3 

0.77 

0.68 

0.49 

0.41 

4 

0.75 

0.42 

0.38 

5 

0.54 

0.49 

6 

0.55 

7 
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evening  emerges  except  for  low  lying  adjacent  layers.  Consider,  for  example,  rn(l,7) 
which  for  successive  evenings  takes  on  the  values  0.44  (12  profiles),  0.08  (39  profiles), 
0.27  (56  profiles),  0.18  (26  profiles),  -0.12  (16  profiles)  and  0.25  (43  profiles),  yielding  a 
weighted  mean  of  0.19.  For  all  192  profiles  rii(l,7)  =  0.37,  that  is,  about  a  factor  of  two 
larger  than  nightly  measurements  would  suggest.  Similarly  for  rn(6,7)  we  have  on  a 
nightly  basis  -0.01,  0.24,  0.60,  0.47,  0.03,  and  0.04,  giving  a  weighted  mean  of  0.30, 
whereas  for  the  entire  sample  we  have  rn(6,7)  =  0.60. 

Statistically  significant  negative  correlations  have  been  obtained  on  some 
evenings.  For  example,  for  2  May,  based  on  39  profiles,  r£^(l,4)  =  -0.31  which  yields  a 
t-statistic^  of  —1.98,  which  is  significant  at  the  0.05  level.  For  4  May  based  on  26 
profiles  we  find  that  rj^(2,7)  =  -0.49  which  yields  a  t-statistic  of  -2.75  which  is 
significant  at  the  0.01  level. 

We  have  also  have  calculated  for  Cn  and  for  log(Cn)  correlation  coefficients  with 
altitude  for  the  entire  data  set,  obtaining  values  of  —0.31  for  Cn  and  -0.50  for  log(Cn)- 

For  each  evening  and  for  the  entire  sample  we  have  calculated  the  slope  and 
intercept  of  the  least  squares  line  for  Cn  and  log(Cn)  with  altitude. 

2.3  Scatter  Plots 

Scatter  plots  have  been  produced  to  illustrate  the  above  discussion.  The  ones 
presented  here  are  all  for  log(Cn)  rather  than  for  Cn  because  including  even  one  outlier 
of  Cn  can  lead  to  a  rather  distorted  plot.  Figure  2  shows  the  graph  of  log(Cn)  for  levels 
1  versus  level  7  for  the  entire  data  set.  This  plot  is  typical  for  a  moderate  amount  of 
correlation.  We  also  show  the  least  squares  line.  Figure  3  shows  a  scatter  plot  for 
log(Cn)  for  level  1  versus  level  2  and  shows  a  high  degree  of  correlation.  In  Figure  4  we 
display  a  typiccil  scatter  diagram  for  a  negative  correlation  coefficient  and  in  Figure  5 
such  a  plot  for  a  near  zero  correlation. 

3.  POSSIBLE  CAUSES  OF  HIGH  CORRELATION 

Various  mechanisms  exist,  other  than  the  scintillometer  itself,  that  could  lead  to 
a  high  degree  of  correlation  of  Cn  and  log(Cn)  values  between  levels.  During  the  course 
of  the  PSU  measurement  program,  atmospheric  conditions  were  not  steady  from 
evening  to  evening,  nor  during  the  course  of  a  single  evening.  Jet  stream  passages  were 

^Brockett,  P.  and  Levine,  A.  (1984),  Statistics  and  Probability  and  Thetr  Applications, 
p.  320,  Saunders  College  Publishing,  1984. 
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LEVEL  1 

Figure  2.  log(Cn  ):  Level  1  versus  Level  7;  30  April-6  May  1986. 
192  profiles.  Shown  for  comparison  is  the  least  squares  line. 
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observed  over  this  time  period  resulting  in  significant  changes  in  the  various  turbulence 
parameters  measured®.  Jet  stream  effects,  for  example,  could  influence  atmospheric- 
conditions  across  most  of  the  scintillometer  measurement  range.  After  looking  at  a 
simple  example  in  Section  3.2  we  estimate  the  possible  effects  of  changing  atmospheric 
conditions  in  two  ways  in  Section  3.3. 

If  instantaneous  localized  measurements  of  Cn  were  available  at  different 
altitudes  we  would  expect  there  to  be  no  correlation.  However,  scintillometer  profiles 
involve  a  good  deal  of  spatial  and  temporal  averaging.  This  effect  will  be  addressed  in 
Section  3.4. 

3.1  Stationarity  and  Independence 

A  stationary  variable®  is  one  which  is  distributed  about  some  definite  mean  with 
a  definite  standard  deviation  and  these  values  do  not  change  with  time.  Two  variables 
are  said  to  be  independent  if  changes  in  one  do  not  influence  the  other.  To  illustrate  the 
effects  of  stationarity  and  independence  (or  the  lack  thereof)  we  construct  the  following 
simple  data  sets. 

Let  X(i,  j)  be  a  random  variable  uniformly  distributed  over  the  interval  [i,  j|. 
X(i,  j)  is  clearly  a  stationary  variable  with  mean  and  standard  deviation  (j  -  i). 
In  Table  9- A,  in  column  one  are  shown  12  values  of  X(l,  3)  and  in  column  two  are 
shown  12  values  of  X(2,  4).  Below  the  table  is  shown  the  correlation  coefficient  for  these 
two  independent  sets.  An  individual  value  of  this  correlation  coefficient  is  inconclusive. 
We  generated  a  hundred  such  tables  and  calculated  the  corresponding  correlation 
coefficients  and  found  an  average  correlation  coefficient  of  —0.01  (with  a  standard 
deviation  of  0.33)  which  is  what  we  would  expect:  the  correlation  between  two 
independent  stationary  variables  should  be  zero.  Tables  9-B  and  9-C  are  self- 
explanatory. 

Consider  Table  10.  It  is  formed  by  merging  Tables  9-A,  B  and  c.  A  high  degree 
of  correlation  is  exhibited.  This  is  not  a  coincidence:  when  100  such  tables  were 
generated,  we  found  that  the  average  correlation  coefficient  was  0.89  with  a  standard 
deviation  of  0.02.  We  can  explain  this  correlation  in  at  least  two  ways: 

®Beland,  R.  (1990),  Geophy.sics  Laboratory,  private  communication. 

®Cox,  D.  R.  and  Lewis,  P.  S.  (1960),  The  Statistical  Analysis  of  Senes  of  Events,  p. 
60,  John  Wiley  &  Sons. 
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Table  9.  Two  sets  of  random  numbers  and  their  correlations,  r. 

A.  The  first  set  is  drawn  from  the  interval  (1,3)  and  the  second  from  (2,4) 

B.  The  first  set  is  drawn  from  the  interval  (3,5)  and  the  second  from  (4,6) 

C.  The  first  set  is  drawn  from  the  interval  (5,7)  and  the  second  fruin  (6,8) 


A 


X(1.3) 

X(2.4) 

1.84 

2.88 

1.32 

2.56 

1.94 

2.31 

2.01 

3.49 

2.36 

3.98 

1.49 

3.74 

1.46 

2.34 

1.61 

2.83 

2.11 

2.99 

2.85 

2.06 

2.85 

2.66 

2.83 

2.66 

r=  -0.14 


B 


X(3.5) 

X{4.6) 

4.76 

4.95 

3.10 

5.38 

3.93 

4.08 

3.57 

4.48 

3.77 

4.15 

3.20 

4.28 

3.93 

msm 

3.94 

3.34 

4.19 

4.59 

4.63 

3.84 

5.83 

3.97 

5.11 

r=  0.11 


C 


X(5,7) 

X(6.8) 

5.92 

7.93 

5.62 

7.95 

6.96 

5.90 

msm 

6.55 

5.33 

■ina 

6.71 

1^1 

6.89 

5.51 

7.58 

6.41 

7.66 

5.92 

6.34 

r  =  -0.08 


(i)  The  two  columns  are  not  independent;  although  each  group  of  twelve  pairs 
were  generated  independently,  these  three  groups  were  not  merged  at  random. 

(ii)  If  we  think  of  each  column  as  a  time  series  with  measurements  made  every 
second,  then  every  12  seconds  the  mean  of  each  set  changes. 

3.2  Possible  Effects  on  r£^(i,  j)^° 

We  think  that  a  good  deal  of  the  observed  correlation  in  Table  8,  especially  for 
non-adjacent  levels,  can  be  explained  in  terms  of  the  above  example.  Atmospheric 
conditions  vaxied  from  one  night  to  the  next  during  the  measurement  program.  Nightly 
means  and  standard  deviations  for  Cn  and  log{Cn)  vary  considerably.  See  Table  11  for 
nightly  values  for  log(Cn).  The  fact  that  these  nightly  values  do  differ  suggests  that  Cn 
is  not  a  stationary  variable  and  that  lumping  all  192  profiles  into  one  data  set  and  then 
calculating  correlation  coefficients  between  levels  is  not  unlike  merging  Tables  9-A,B 
and  C  to  form  Table  10.  In  order  to  see  the  type  of  correlation  that  may  be  due  to  this 
phenomenon  we  take  two  approaches. 


^°Because  rn(i,  j)  and  1^(1,  j)  are  comparable  in  magnitude,  we  concentrate  on  the 


latter. 
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Table  10. 

The  result  of  merging 
Tables  9A,  B  and  C 


3.2.1  Treat  Cn  as  a  ‘‘Nightly”  Random  Variable 

Suppose  we  assume  that  on  any  given  night  Cn  has  a  uniform  distribution  whose 
mean  and  standard  deviation  at  each  level  are  as  given  in  Table  11.  This  assumption 
will  underestimate  the  amount  of  correlation  from  level  to  level  in  that  Cn  is  most 
likely  log  normally  distributed^^.  Furthermore,  during  the  course  of  an  evening, 
atmospheric  conditions  may  have  changed.  We  now  take  the  nightly  mean  and  standard 
deviation  for  log(Cn)  and  generate  n  random  variables  (where  n  is  the  number  of 
profiles  for  the  night  in  question)  centered  on  this  mean  in  a  range  of  width  2^13  times 
the  nightly  standard  deviation^^.  We  then  combine  results  for  each  together  as  in  the 
above  example  to  get  192  profiles.  We  then  calculate  correlation  coefficients  between 
levels  using  this  data.  We  then  repeated  this  procedure  100  times.  Table  12  shows  the 
resulting  average  correlation  coefficients.  For  non-adjacent  levels  especially,  the 
resulting  correlations  are  comparable  to  those  of  Table  8. 

Table  11.  Nightly  means  and  standard  deviations  for  log(Cn)  by  level.  The  units  of 
_2 

Cn  are  m  ^  and  n  is  the  number  of  profiles. 


DATE 

■ai 

1 

2 

3 

4 

5 

6 

7 

4/30 

■ 

-16.123 

0.538 

-16.749 

0.316 

-16.935 

0.180 

-16.948 

0.192 

-17.187 

0.191 

-17.117 

0.196 

5/2 

39 

-16.107 

0.156 

-16.512 

0.135 

-16.944 

0.291 

-16.934 

0.197 

-17.050 

0.163 

-17.023 

0.133 

5/3 

56 

-15.800 

0.426 

-16.709 

0.306 

-17.117 

0.339 

-17.077 

0.232 

-17.138 

0.223 

-17.198 

0.197 

5/4 

26 

-16.599 

0.130 

-17.130 

0.266 

-17.548 

0.331 

-17.639 

0.412 

mm 

-17.271 

0.279 

5/5 

16 

-16.790 

0.184 

-17.350 

0.221 

-17.853 

0.446 

m 

-17.766 

0.407 

-17.566 

0.146 

5/6 

43 

-16.661 

0.601 

-16.796 

0.488 

-17.065 

0.427 

■a 

-17.473 

0.191 

^^Loos,  G.  and  Hogge,  C.  (1979),  Turbulence  of  the  Upper  Atmosj)licre  and 
Isoplanatism,  Appl.  Opt.,  18:  2654-2661. 

^^One  can  show  that  a  uniform  distribution  centered  on  x  with  width  2  ■js  .s  will  have 
a  standard  deviation  of  s. 
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3.2.2  Recalculate  Table  8  with  Nightly  Means  Subtracted 

A  second  approach  in  dealing  with  the  non-stationary  nature  of  Cn  is  to  sulrtract 

out  the  nightly  means  before  calculating  correlation  coefficients.  The  result  of  doing  this 

is  shown  in  Table  13.  This  table  should  again  be  compared  to  Table  8.  Correlation  has 

been  significantly  reduced.  For  example,  7)  has  dropped  from  0.52  down  to  0.20. 

The  amount  of  variation  in  log(Cn)  at  level  7  that  can  be  inferred  from  variation  at 

o 

level  1  has  dropped  from  27  percent  down  to  rj^(l,  7)  =  4  percent.  We  further  note  that 
r^(l,  4)  =  8  percent  and  r£(4,  7)  =  2  percent  which  would  indicate  only  a  slight 
interdependence  between  levels  1,  4,  and  7.  This  residual  correlation  may  be  due  to  the 
scintillometer  but  could  also  be  due  to  non-stationarity  of  Cn  on  a  given  evening. 

3.3  Altitude  Dependence  of  Cn 

Various  models  for  Cn  have  been  proposed.  In  most  of  these  altitude  is  directly 
or  indirectly  the  dominant  parameter.  Any  model  of  the  form  Cn(z)  =  Cn(zo)  f(z) 
w’here  f(z)  is  any  function  of  altitude  and  Cn(zo)  is  some  reference  value  would  predict  a 
high  degree  of  correlation,  but  not  necessarily  linear  correlation.  Some  such  models 
include  that  of  Kairnal  as  modified  by  Walters  and  Kunkel,^^  that  of  Battles,  Murphy 
and  Noonan^'*  and  that  of  Kukharets  and  Tsang,^^  each  of  which  would  yield  a  high 
correlation  of  Cn  from  level  to  level.  We  are  not  saying  that  these  models  do  apj)ly  to 
the  altitude  range  covered  by  this  data  but  are  suggesting  that  the  high  degree  of 
correlation  observed  may  be  due  to  the  fact  that,  although  very  different  physical 
mechanisms  are  producing  turbulence  at  the  different  levels  sampled,  when  sufficient 
data  are  used  we  expect  that  on  the  whole  Cn  will  be  a  function  of  altitude.  This  would 
explain  the  difference  between  nightly  results  versus  that  observed  over  the  entire 
measurement  program.  From  night  to  night  differing  atmospheric  conditions  exist. 
When  these  are  averaged  over  all  six  nights  Cn  becomes  correlated  due  to  the  fact  that 
it  is  overall  a  decreasing  function  of  altitude. 

The  correlation  coefficients  calculated  for  Cn  and  log(Cn)  with  altitude  are  a 
further  indication  that  altitude,  not  instrumentation,  is  affecting  the  values  of  ri^(i,j) 

^^Walters,  D.  L.  and  Kunkel,  K.  E.  (1981),  Atmospheric  Transfer  Function  for  Desert 
and  Mountain  Locations:  The  Atmospheric  Effects  on  Cq,  J.  Opt.  Soc.  Am.,  71;  397-405. 

^‘’Battles,  F.  P.,  Murphy,  E.  A.  and  Noonan,  J.  P.  (1988),  The  Contribution  of 
.Atmospheric  Density  to  the  Drop-off  Rate  of  Cn,  Physica  Sertpta,  37:  151-153. 

^^Kukharets,  V.  P.  and  Tsang,  L.  R.  (1980),  Structure  Parameter  of  the  Refractive 
Index  in  the  Atmospheric  Boundary  Layer,  IZV  Acad  Sci  USSR,  Alrno.s.  Oceanic  Pliys.,  16; 
397-405. 
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Table  12.  Correlation  coefficients  obtained  by  assuming  that  log(Cn)  is  randomly 
distributed  with  a  mean  and  standard  deviation  as  given  in  Table  11.  These 
results  are  based  on  the  average  of  100  simulations  and  each  has  a  standard 
deviation  of  about  0.06. 


1 

2 

3 

4 

5 

6 

7 

1 

0.45 

0.50 

0.39 

0.47 

0.46 

0.41 

2 

0.50 

0.37 

0.43 

0.37 

0.26 

3 

0.44 

0.48 

0.40 

0.34 

4 

0.40 

0.33 

0.29 

5 

0.46 

0.38 

6 

0.39 

7 

Table  13.  Correlation  coefficients  for  log(Cn)  obtained  from  the  scintillometer 
data  after  the  nightly  means  given  in  Table  11  have  been  subtracted  out. 


1 

2 

3 

4 

5 

6 

7 

1 

0.79 

0.62 

0.28 

0.29 

0.20 

2 

0.83 

0.44 

0.30 

0.14 

3 

0.64 

0.38 

0.20 

0.14 

4 

0.61 

0.16 

0.15 

5 

0.20 

0.19 

6 

0.24 

7 

22 


and  rj^(i,j).  Should  a  model  of  the  above  suggested  form  hold  we  would  expect  that  a 
least  squares  line  for  log(Cn)  from  level  to  level  would  have  a  slope  of  near  unity.  Our 
calculations  for  adjacent  levels  yield  slopes  of  order  0.8. 

If  the  observed  correlations  were  primarily  instrument  caused,  we  would  expect  a 
similar  pattern  on  any  given  evening  similar  to  that  exhibited  overall:  statistically 
significant  positive  correlation  that  decreases  with  level  separation.  Weighting  function 
overlap  can  only  result  in  positive  correlation,  yet  significant  negative  correlation  has 
been  observed  on  several  occasions. 

4.  CONCLUSIONS 

Because  rL(i,  j)  can  be  interpreted  as  giving  the  percent  variation  in  log(Cn)  at 
level  i  that  is  explained  by  least  squares  regression  on  log(Cn)  at  level  j,  we  present 
Tables  14-16  which  give  the  values  of  r£(i,  j)  for  the  192  profiles  as  measured  (the 
square  of  Table  8),  the  values  obtained  by  assuming  that  log(Cn)  is  uniformly 
distributed  about  the  nightly  mean  (the  square  of  Table  12)  and  the  values  of  1-^(1,  j) 
obtained  by  subtracting  out  the  nightly  means  (the  square  of  Table  13.) 

Suppose  that  we  somewhat  arbitrarily  define  “substantial  correlation”  to 
correspond  to  r£(i,  j)  >  0.20  and  “insubstantial  correlation”  to  correspond  to 

9  .  . 

rL(b  j)  <  0.10.  Then  in  Table  14,  we  see  that  there  is  substantial  correlation  between 
practically  all  pairs  of  levels.  However,  when  we  consider  Table  16  we  see  that 
correlation  between  almost  all  non-adjacent  levels  is  insubstantial,  the  exceptions  being 
rL(l,  3)  =  0.38,  rL(2,  4)  =  0.19  and  rL(3,  5)  =  0.14.  There  is  insubstantial  correlation 
shown  for  all  of  the  r£(i,  j)  for  i  <  j  -f  1  when  the  nightly  means  have  been  removed. 

In  Table  15  we  see  that  the  non-stationarity  of  Cn  may  have  introduced 
substantial  correlation  between  some  levels,  especially  when  it  is  remembered  that  this 
table  was  generated  using  assumptions  that  underestimate  the  values  of  rL(i,  j)  as  noted 
in  Section  3.2.1. 

We  conclude  that  values  obtained  for  Cn  at  levels  1,  4,  and  7  are  relatively 
independent  of  one  another.  This  also  appears  to  be  true  for  all  but  low  lying  adjacent 
or  next-adjacent  levels. 


23 


Table  14.  The  square  of  the  correlation  coefficients  between  levels  for  log(Cn). 
April  30-May  6,  1986  (192  profiles).  Compare  with  Table  8. 


Table  15.  Square  of  the  correlation  coefficients  obtained  by  assuming  that 
log(Cn)  is  randomly  distributed  with  a  mean  and  standard  deviation  as  given  by 
Table  11.  These  results  are  based  on  the  average  of  100  simulations.  Compare 
with  Table  12. 
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